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Dedicated research into the design and construction of a large scale Quantum Information Pro- 
cessing (QIP) system is a complicated task. The design of an experimentally feasible quantum 
processor must draw upon results in multiple fields; from experimental efforts in system control and 
fabrication through to far more abstract areas such as quantum algorithms and error correction. 
Recently, the adaptation of topological coding models to physical systems in optics has illustrated 
a possible long term pathway to truly large scale QIP. As the topological model has well defined 
protocols for Quantum Error Correction (QEC) built in as part of its construction, a more grounded 
analysis of the classical processing requirements is possible. In this paper we analyze the require- 
ments for a classical processing system, designed specifically for the topological cluster state model. 
We demonstrate that via extensive parallelization, the construction of a classical "front-end" system 
capable of processing error correction data for a large topological computer is possible today. 

PACS numbers: 03.67.Lx, 07.05.Wr 



I. INTRODUCTION 

The design and construction of a feasible large scale 
quantum computing system has been a highly sought af- 
ter and long term goal of quantum information science 
ever since the first physical system proposals were made 
in the mid 1990's [IlEliaiaElElIIlEllSlIIO]. While 
experimental advances in quantum computing have been 
pronounced [III [11 111 [11 [13 [11 [171 [^1 [E] we are 
yet at the stage where we can faithfully claim a multi- 
million qubit device is just around the corner. 

Nevertheless, in order for experimental progress to be 
made, the fundamental theoretical building blocks for a 
large scale computer need to be firmly in place. This the- 
oretical development is not restricted to the discovery of 
new protocols for computation, algorithms or error cor- 
rection but it also includes the architectural engineering 
of future computers. 

While there has been steady progress over the past 15 
years on designing novel and (more importantly) experi- 
mentally feasible large scale processor architectures, the 
complication of implementing appropriate and efficient 
error correction procedures and designing systems which 
can trivially be scaled to the level of a "programmable" , 
multi-task computer is still a daunting and often ne- 
glected area of research. 

Recently, the introduction of theoretical ideas 
such as topological cluster state quantum computing 
(TCQC) [21 [m [22] and the single photon optical ar- 
chitecture [21 |2l] gives us an idea of what a truly large 
scale device may possibly look like. The modular design 
of the cluster preparation network, and the measurement 
based nature of the computational model, gives this de- 
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sign something that other architectures arguably lack, a 
strictly modular scaling of the entire computer. 

While the design introduced in Refs. [23' ^ is not 
necessarily the optimal way to construct a large scale 
quantum computer, it does contain several key elements, 
easing the conceptual design of a large scale computer. 
For example: 

1. Utilizing a computational model that is fundamen- 
tally constructed from error correction, rather than 
the implementation of codes on top of an otherwise 
independent computational model. 

2. Having a modular construction to the computer. 
The fundamental quantum component (the pho- 
tonic chip [2l [26] ) , is a comparatively simple quan- 
tum device. Scaling the computer arbitrarily re- 
quires the addition of more chips in a regular and 
known manner. 

3. Employing a computational model exhibiting high 
fault-tolerance thresholds [20], which relieves the 
pressure on experimental fabrication and control. 

4. Utilizing a measurement based model for compu- 
tation [27 . By employing a measurement based 
computational model, the quantum component of 
the computer is a simple state preparation network. 
Therefore, programming such a device is a problem 
of classical software, not of hardware. 

These properties, as well as others, allowed us to con- 
sider the structure of an extremely large mainframe-type 
device. In Ref. [25 , the quantum analogue of high per- 
formance computing was examined. The conceptual scal- 
ability of the optical architecture allowed us to examine 
the operating conditions, physical structure and resource 
costs of a computer employing extensive topological error 
correction to the level of 2.5 million logical qubits. 
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In addition to examining the implementation of a large 
scale TCQC, the nature of the topological model also 
allows for a more concrete discussion on an often ne- 
glected, but important aspect of quantum information 
processing, namely what are the classical computational 
requirements of a large scale device? In this paper we 
attempt to answer this question. 

There have been several broad investigations into the 
classical structure, design and operation of a large scale 
quantum computer [28^, but investigation into this 
topic is difficult. The primary obstacle in analyzing clas- 
sical requirements is that the quantum architecture gen- 
erally has to be specified. Since all classical processing is 
ultimately dependent on both the computational model 
employed at the quantum level and more importantly the 
error correction protocols utilized, a detailed analysis of 
the classical front end must wait for the design of the 
quantum processor. 



In this paper we specifically analyze the classical front 
end requirements to perform active error correction on a 
3D topological cluster lattice prepared by the photonic 
chip network. This analysis will be restricted to the clas- 
sical system required to implement the underlying error 
correction procedures in the topological model, without 
the execution of an active quantum algorithm. 

Although we present this analysis in the context of 
the optical network presented in Ref. [23 , it should be 
stressed that this analysis is still highly relevant for any 
physical architecture employing the 2D or 3D topologi- 
cal model [201 El EQI ISTJ |32 . Our analysis demonstrates 
that with several optimizations of the classical process- 
ing and the ability to significantly parallelize classical 
error correction processing, the classical computational 
requirements for large scale TCQC are indeed within the 
capabilities of today processing technology. 

Section |n| very briefiy reviews the nature of the topo- 
logical cluster model in order to fully specify what is re- 
quired of the classical processing. Section [ill] reviews the 
fiowing nature of the preparation network, and how this 
relates to the optical measurement layer, and the first 
level of classical processing. In Section IV we overview 



the basic requirements of the classical network and how 
target processing rates are related to the clock cycle of 
the quantum network. In Section [V| we introduce classi- 
cal benchmarking data for the minimum weight match- 
ing algorithm utilized for error correction and illustrate, 
given this data, how error correction processing can be 
parallelized over the entire computer. We conclude by il- 
lustrating how the parallelization of the classical process- 
ing allows, in principal, for the construction of a classical 
error correcting front end for large scale TCQC with clas- 
sical processing technology available today. 



II. TOPOLOGICAL ERROR CORRECTION IN 
THE OPTICAL ARCHITECTURE 

TCQC was first introduced by Raussendorf, Harring- 
ton and Goyal in 2007 [201 EI]. This model incorporates 
the ideas stemming from topological quantum comput- 
ing introduced by Kitaev [33] and cluster state compu- 
tation [27 , leading to a very attractive computational 
model incorporating error correction by construction and 
exhibiting a high threshold error rate. 

As with any measurement based computational model, 
computation proceeds via the initial construction of a 
highly entangled multi-qubit state. Fig. [l] illustrates 
the structure of the cluster. Each node in the clus- 
ter represents a physical qubit, initially prepared in the 
1+) = (|0) + |1))/a/2 state, and each edge represents a 
controlled-cr^ entangling gate between qubits. This is the 
fundamental unit cell of the cluster, which repeats in all 
three dimensions. 

Computation under this model is achieved via the con- 
sumption of the cluster along one of the three spatial 
dimensions [22 (simulated time). Logical qubits are de- 
fined via the creation of "holes" or "defects" within the 
global lattice and multi-qubit operations are achieved via 
braiding (movement of these defects around one another) 
as the cluster is consumed along the direction of simu- 
lated time. The specific details for computation under 
this model are not important for this discussion and we 
encourage the reader to refer to Refs. |20l [22] for fur- 
ther details. For this analysis, the effect of errors on a 
topological lattice is the important factor. 



A. Error effects 

Quantum errors in this model manifest in a very spe- 
cific way. In Fig.[l] illustrating the unit cell of the cluster, 
we have illustrated six face qubits shown in red. If no er- 
rors are present in the system, measuring each of these 
six qubits in the |+) or |— ) states {ax basis) will return 
an even parity result. If we denote the classical result of 
these six measurements as Si G {0, 1}, i = 1, .., 6, then 



{Si + 52 + 53 + S4 + 55 + Sq) mod 2 = 0. 



(1) 



This result is a consequence of the quantum correlations 
established in the preparation of the cluster. If the cluster 
is prepared perfectly, these six qubits are placed into a 
state that is a +1 eigenstate of the operator 



K = al 



(2) 



where is the Pauli bit-fiip operator. The measurement 
of each of these six qubits in the basis will produce 
random results for each individual qubit. However the 
eigenvalue condition of this correlation operator guaran- 
tees the classical measurement results satisfy Eq.[l]in the 
absence of errors. 



FIG. 1: Unit cell of the 3D cluster, which extends arbitrarily 
in all three dimensions. In the absence of errors, the entan- 
glement generated within the lattice sets up correlation con- 
ditions when certain qubits are measured. For each unit cell 
of the cluster, quantum correlations guarantee that the mea- 
surement results of the six face qubits (shown above) return 
an even parity result [Eq. [l] . 



The remaining qubits in each unit cell are also mea- 
sured in the Gx basis, but their results are associated 
with the parity of cells within the dual lattice [Fig. 2]. 
This property of the cluster is not important for this dis- 
cussion. What is important is that when no quantum 
algorithm is being implemented, every qubit is measured 
in the cf^ basis and is used to calculate the parity of their 
respective cells. 

Due to the structure of this 3D cluster state, all er- 
ror channels can effectively be mapped into phase errors 
(Pauli-cr^ operations applied to phyical qubits in the lat- 
tice which takes a state, a|0) + Q^|0) — I3\l)) or 
physical qubit loss [22j. These two distinct channels are 
processed slightly differently. 



B. Error channel one: Phase errors 

We first consider the effect of phase errors, which act 
to flip the parity of a cell of the cluster. As the Pauli 
operators (7^ and cr^ anti-commute, (JzC^x = —(^x^zi if a 
phase error occurs to one of these six qubits the correla- 
tion condition of a cell will flip from a +1 eigenstate of 
the operator to a -1 eigenstate. If the correlation con- 
dition flips to a —1 eigenstate of the classical result 
of the six individual measurements also flips to an odd 
parity condition, i.e. 



{Si + 52 + S3 + 54 + 55 + 56) mod 2 = 1. 



(3) 



Errors on qubits on the boundary of each unit cell there- 
fore flip the parity of the measurement result from even 



FIG. 2: (From Ref. [22]) The regular structure of the 3D clus- 
ter results in a primal and dual lattice structure. A set of eight 
unit cells arranged in a cube results in a complete unit cell 
present at the intersection of these eight cells. Hence there are 
two self-similar lattices (offset by half a unit cell diagonally) 
known as the primal and dual lattice. These two structures 
are extremely important for computation under the topologi- 
cal model [20]. However, in the context of this discussion, the 
measurement results of the additional nine qubts in Fig.[l]are 
associated with the parities of bordering dual cells. 



to odd. Note that the change of parity of any individual 
cell gives us absolutely no information regarding which 
one of the six qubits experienced the physical error. 

The second important aspect of this structure is that 
any given qubit lies on the boundary of two cells in the 
cluster. If a given qubit experiences a phase error it will 
flip the parity result of the two adjacent cells. This al- 
lows us to detect which of the six qubits of a given cell 
experienced an error. If a single cell flips parity, we then 
examine the parity result of the six adjacent cells. As- 
suming that only one error has occured, only one of these 
six adjacent cells will have also flipped parity allowing us 
to uniquely identify the erred qubit. 

If we now consider more than one error within the lat- 
tice, we no longer identify the location of individual errors 
but instead identify error chains. Fig. |3] from Ref. [22] 
illustrates. Here we have a 3D cluster consisting of a 
cube of 4^ cells with three error chains. The first chain is 
a single error which flips the parity of the two adjacent 
cells, the other two chains illustrate the effect of multiple 
errors. As the parity conditions are cyclical (mod 2), if 
two errors occur on the boundaries of a given cell the 
parity result will not change. Instead, the two cells at 
the endpoints of these error chains are the cells which 
flip parity. 

Hence, in the TCQC model, it is not the locations of 
individual errors which are important but the endpoints 
of error chains. In fact, the symmetries of the cluster do 
not require us to identify the physical error chain corre- 



FIG. 3: (Taken from Ref. [22]) Illustration of error effects in the 3D cluster. Here we show a volume 1/ = 4^ of cluster cells 
and the effect of three error chains. As the parity conditions of Eq. [ijare cyclical (mod 2), the calculation of cell parities only 
reveals information regarding the endpoints of error chains. Shown above are three examples, one error flips the parity of the 
two cells adjacent to the erred qubit while longer chains only flip the parity of the end point cells. The goal of error correction 
is to faithfully identify these chains given the end point parity data. 



spending to the detected endpoints. Once the endpoints 
of the chain are correctly identified all path of correc- 
tion operators (Pauli operators which are applied to re- 
verse detected errors) which connect the two endpoints 
are equivalent [20 . Hence, the goal of error correction in 
this model is to correctly "pair up" a set of odd parity 
cells such that the appropriate correction operators can 
be applied. 

Undetectable errors in this model occur when chains 
become so long that they actually connect two bound- 
aries of the lattice. If a physical error chain completely 
spans the lattice from one boundary to another then each 
individual cell experiences two physical errors and every 
cell remains in an even parity state. If the 3D lattice is 
not used for computation, these error chains are actually 
invariants of the cluster and hence have no effect. Once 
computation begins, information is stored by deliberately 
creating holes (or defects) in this lattice. These defects 
act as artificial boundaries and consequently error chains 
connecting defects to any other boundary (either other 
defects or the boundary of the lattice) are undetectable 
and cause logical errors on stored information. 

From the standpoint of this investigation we are only 
concerned with performing active error correction on a 
defect free lattice. We will not be introducing informa- 
tion qubits into the cluster. Instead we will be examining 



the classical resources required to detect and correctly 
identify error chains in an otherwise perfectly prepared 
lattice. 

This type of analysis is justified as information qubits 
are essentially regions of the 3D cluster that have sim- 
ply been removed from the global lattice. Analyzing the 
classical requirements for the complete, defect free lat- 
tice therefore represents the maximum amount of classi- 
cal data that needs to be processed for correction. 



C. Error channel two: Qubit loss 

The second major error channel is qubit loss. As we 
are motivated by the optical architecture introduced in 
Ref. [23], loss represents a significant error channel. Un- 
like other computation models for quantum information, 
the TCQC model can correct loss without additional ma- 
chinery. Without going into the details, loss events can 
be modeled by tracing out the qubit that is lost from the 
system and replacing the qubit with the vacuum. Tracing 
out the lost qubit is equivalent to measuring the qubit in 
the |0) or |1) state (a basis measurement) with the 
result unknown. In principal, this type of error can be 
modeled as a standard channel, causing the parity of the 
respective cells to flip with a probability of 50%. How- 
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ever, since the qubit is no longer present, we can utilize 
the vacuum measurement to uniquely identify these error 
events. 

Illustrated in Fig. [4] is the structure of a unit cell when 
one qubit is essentially measured out via loss. In this 
case, the boundary of a cell increases around the lost 
qubit. Instead of the parity conditions being associated 
with the six face qubits of a given cell, it extends to be 
the combined parity of the ten measurements indicated. 
As the loss event is detected via no "clicks" from the 
detector array, this result is corrected by now taking the 
parity of this larger structure and proceeding as before. 
Provided no other errors have occured, the parity of this 
larger boundary will be even, and any additional qubit 
errors will link this extended cell to a second end point 
with odd parity. Recent results, obtained in the context 
of the surface code [STJ [32], have demonstrated a high 
tolerance to heralded loss events [34j. 




FIG. 4: The effect of qubit loss on the parity conditions of 
the cluster. When a qubit is lost, it is essentially removed 
from the cluster. The parity condition of Eq.[l]is extended to 
the boundary surrounding the loss event. As a lost qubit is 
a heralded error (i.e. can be detected separately), the parity 
calculation can be modified to encapsulate this larger volume 
if a loss event is detected. 



III. THE OPTICAL COMPUTER, A 
"FLOWING" NETWORK 

As introduced in Ref. [23 , optical TCQC can be per- 
formed by making use of a preparation network of pho- 
tonic chips [26 . This network receives a stream of sin- 
gle photons (from a variety of appropriate sources) and 
ejects a fully connected topological cluster. As compu- 
tation proceeds via the successive consumption of lattice 
faces along the direction of simulated time, the prepa- 
ration network is designed to continuously prepare the 



lattice along this third dimension at a rate equal to the 
rate of consumption by the detector layer. Fig. |5] illus- 
trates the basic design. 

Photons are continuously injected into the rear of 
the preparation network, ideally from appropriate on- 
demand sources. Each photon passes through a network 
of four photonic chips, which act to link them together 
into the appropriate 3D array. Each photonic chip op- 
erates on a fundamental clock cycle, T, and each chip 
in the network operates in a well-defined manner, inde- 
pendent of the total size of the network [23]. In total, a 
single photon entering the network at t = will exit at 
t = 4T, after which it can be measured by the detector 
banks. 

Each photonic chip acts to entangle a group of five 
photons into an appropriate state such that the parity 
condition for each cell is satisfied. After each group of 
five photons passes through an individual chip, a sin- 
gle atomic system contained within each chip is mea- 
sured and reinitialized, projecting the relevant group of 
5-photons into an entangled state. The result of this 
atomic measurement is fed forward to the classical pro- 
cessing layer in order to define a set of initial correlation 
conditions. The cluster is defined such that Eq{l]is satis- 
fied for all cells. However, the preparation network does 
not automatically produce these correlation conditions. 
Depending on the measurement results of the atomic sys- 
tems contained within each photonic chip, approximately 
50% of cells within the lattice will be prepared with an 
initial parity condition that is odd. This can, in principal, 
be corrected to be even for all cells by applying selective 
single qubit rotations dependent on the atomic readout 
results, but this is unnecessary. The initial parity results 
from the preparation network are simply recorded, end- 
points of error chains are then identified with cells that 
have changed parity from this initial state. 

As one dimension of the topological lattice is identi- 
fied as simulated time, the total 2D cross section defines 
the actual size of the quantum computer. Defects, re- 
gions of the cluster measured in the basis, are used 
to define logical qubits and are kept well separated to 
ensure fault tolerance. The 2D cross section is then con- 
tinually teleported, via measurement, to the next succes- 
sive layer along the direction of simulated time allowing 
an algorithm to be implemented (in a similar manner to 
standard cluster state computation [27 ). 

In Fig. |6] we illustrate the structure of the detection 
system. For the sake of simplicity, we are assuming that 
the basis states for the qubits are photon polarization, 
\H) = |0) and |V) = hence the detection system con- 
sists of a polarizing beam splitter (PBS) and two single 
photon detectors. A given unit cell flows through a set 
of nine optical lines such that the relevant parity is given 

by, 
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FIG. 5: General architectural model for a "flowing" optical computer. Single photon sources inject photonic qubits into a 
preparation network that deterministically links up a 3D photonic cluster, layer by layer. Immediately after the preparation 
network, an array of single photon detectors measures each photon to perform computation. As photons are continually linked 
to the rear of the 3D cluster as the front is consumed, an arbitrarily deep 3D cluster can be prepared and consumed with finite 
space. 



r 



(4) 



where sj^^ is the detection result (1, 0) of detector (i, j) at 
time T. This result defines the parity of the cell (i, j, T) 
in the lattice. Loss events would result in neither detector 
firing, at which point the calculation of Eq. [4] would be 
redefined based on the loss event to calculate the parity 
of the boundary around this lost qubit. 

The results of all the detection events are fed directly 
from the detectors into the first classical processing layer. 
This layer calculates Eq. [4] and passes the result forward 
to the subsequent processing layer if it differs from the 
initial parity. This general structure extends across the 
entire 2D cross section of the lattice with parities re- 
peatedly calculated for each unit cell that flows into the 
detection system. 



IV. CLASSICAL PROCESSING 
REQUIREMENTS 

Illustrated in Fig. [7| is the layer structure of the clas- 
sical processing for topological error correction. In total 
there are four stages, parity calculation, tree creation. 



minimum weight matching, and the quantum controller. 
As this analysis is only considering the classical require- 
ments for base level error correction processing, we will 
not discuss the structure for the quantum controller. 
This top-level processor will be responsible for the ap- 
plication of active quantum algorithms on the topologi- 
cal lattice, given inputs from both a quantum algorithm 
compiler, and error data outputted from the correction 
processors. In the future, we will be examining the re- 
quirements for this top level processor, but for now we 
omit details regarding its structure and operation. 

The error correction in this model requires identifying 
all cells which have changed parity and reliably identi- 
fying pairs of parity changes associated with a physical 
error chain. In order for error correction to be effective, 
we assume a standard, stochastic, qubit error model. As 
all standard error channels (all errors except qubit loss) 
can be mapped under the topological model to phase er- 
rors [20l[22|, we can, without loss of generality, assume 
that each qubit experiences a phase error, cr^ = Z, with 
probability p. Therefore, a stochastic hierarchy exists 
with respect to longer error chains. 
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FIG. 6: basic structure of the detection layer in the optical network. A given unit cell of the cluster can be associated with 9 
optical waveguides containing temporally separated photonic qubits. For each set of 9 detectors (consisting of a polarizing beam 
splitter and two single photon detectors for polarization encoding), the central detector is associated with the cross-sectional 
co-ordinate address, (i, j), for the cell. The temporal co-ordinate, T, is associated with the current clock cycle of the quantum 
preparation network. For each unit cell, in the absence of photon loss, the measurement results for the 6 face qubits are sent 
to the first classical processing layer that calculates Eq. [4]for each 3D co-ordinate, T). If the parity result differs between 
co-ordinates {i,j,T — 2) and {i,j,T), this information is sent to the next classical processing layer. 



As the probability of a single phase error is given by 
p < 1, the probability of an d error chain is 0{p^) <C 1. 
Therefore, for a given set of parity results, the most likely 
set of physical errors producing the classical measure- 
ment pattern is the set of end point pairings where the 
total length of all connections is minimized. 

Classical results stemming from the detection layer are 
used to calculate the parity for all unit cells for some total 
volume. The co-ordinates of all cells which have experi- 
enced a parity flip are stored in a classical data structure. 
Minimum weight matching algorithms [35l [36[ |37l |38^ are 
then used to calculate likely error chains corresponding 
to the detected endpoints. Once chains are calculated, 
the Pauli frame (the current eigenvalue condition for all 
cells relative to their initial state) of the computation, 
within the quantum controller, is updated with the new 
error information. 

The frequency of error correction in the lattice is dic- 
tated by the application of the non-Clifford gate 



This gate is required to generate a universal set of op- 
erations [39|, and in the topological model is achieved 



via the injection of multiple, low fidelity, singular an- 
cilla, magic-state distillation [40^, "4T and the application 
of teleportation protocols with information qubits. In 
order to successfully apply logical T gates, error infor- 
mation must be obtained for the qubit undergoing the 
teleported T gate prior to application of the correction. 
Fig. [8] illustrates the teleportation protocol to implement 
an RziO) rotation. If a logical X error exists on the state 
prior to teleportation, the condition 

RMx\^) = xR,{-em (6) 

implies that this error must be known before teleporting 
the rotation Rz{0). If the error is detected after telepor- 
tation, the conjugate rotation Rz{—0) will actually be 
applied. Therefore, the classical processing of the mini- 
mum weight matching algorithm will have to occur at a 
comparable rate to the logical gate rate of the prepara- 
tion network to ensure up-to-date error information for 
all logical qubits is available when teleported gates are 
applied. 

As detailed in Ref. [42], the clock cycle of the prepa- 
ration network can vary from nanoseconds to microsec- 
onds, depending on the system utilized to construct the 
photonic chip. Hence our goal in this investigation is to 
determine, for a given failure probability of the quantum 
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FIG. 7: General cross-sectional processing structure for the error correction procedures in the optical TCQC. The quantum 
network consists of a bank of single photon sources, photonic chips and detectors. The preparation network continually prepares 
the cluster along the dimension of simulated time. The detector array will be outputting classical results on the same time frame 
as the quantum clock cycle. The classical processing layer is required to determine likely physical chains of errors producing 
the classical parity results detected. The classical layer consists of 4 stages, parity check, co-ordinate tuple tree creation, 
minimum weight matching, and the quantum controller. 
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FIG. 8: Standard teleportation circuit required to perform 
the rotation Rz{0) on an arbitrary quantum state via the 
preparation of an appropriate ancilla state. The presence of 
a bit-flip = X errors on the qubit affects the gate. Error 
information during quantum processing must be up to date 
upon applying these types of gates in order to ensure rotations 
are applied in the correct direction. 



V. LAYERS TWO AND THREE: MINIMUM 
WEIGHT MATCHING 

Calculating the minimum weight matching of the clas- 
sical parity data stemming from the first layer of the 
classical processing network is the essential requirement 
for error correction in the topological model. The parity 
processing layer is designed to simply output co-ordinate 
tuples for all parity changes in the lattice to this next 
layer. The relevant question is, can the minimum weight 
processing of this data be performed over a large volume 
of the cluster in a comparable time frame to the quantum 
preparation network? 



component of the computer, how quickly the network be 
operated such that all classical processing can be per- 
formed utilizing today's technology. 



A. Minimum Weight Matching benchmarking 

Classical algorithms for determining the minimum 
weight matching of a connected graph are well known 
with a run-time polynomial in the total number of 
nodes 071 |38]. Such algorithms are derived from the 
original Edmonds' Blossom algorithm [36] and for our 
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benchmarking tests we have used Blossom V [38]. How- 
ever, due to the nature of our problem, there are several 
adaptations that can be made to optimize the algorithm 
further. 

Typical minimum weight matching algorithms accept 
a list of N co-ordinates, N — even, and of weighted edges 
(in this case, lattice separation between nodes) such that 
the corresponding graph is completely connected. The 
output is then a list of edges such that every node is 
touched by exactly one edge and the total weight of the 
edges is a minimum. For the purposes of TCQC we can 
optimize by considering the specifics of the problem. 

Due to the stochastic hierarchy of errors in the qubit 
model, and the assumption that the operational error 
rate of the computer is low, p <C 1, the most likely pat- 
terns of errors are simply sets of sparse single errors caus- 
ing two adjacent cells to fiip parity. The probability of 
obtaining longer and longer errors chains is increasingly 
unlikely. Therefore, for a given volume of classical par- 




FIG. 9: Likely structure of the error data for a large volume 
of cluster cells. As we are assuming the operational error rate 
of the quantum computer is low, p ^ 1, long error chains 
become increasingly unlikely. Hence, cell parity flips will most 
likely be sparse and, generally, pairs of parity flips will be 
isolated. This property of the computational model allows a 
certain amount of optimization of the classical requirements 
for minimum weight matching. 

ity results, erred cells will tend to be clustered into small 
groups, as illustrated in Fig. [9] Additionally, we can con- 
sider the computational lattice structure. Fig. [Toj taken 
from Ref. [23 , illustrates the 2D cross-section of the cell 
structure once logical qubits are defined. In this exam- 
ple, defects are separated from each other (and from the 
edge of the lattice) by a total of 16 cells. Logical errors 
occur when error chains connect two boundaries in the 
cluster. If an error chain spans more than 8 cells, the 
correction inferred from the endpoints is likely to be in- 




20 

FIG. 10: (Taken from Ref. [23]). Cross section of a large 
topological lattice. Qubits are defined within a cluster region 
of approximately 40 x 20 cells. The actual qubit information 
is stored in pairs of defect regions (artificially created holes 
in the lattice). As undetectable error chains occur when er- 
rors connect two boundary regions, topological protection is 
achieved by keeping defect qubits well separated from bound- 
aries and from each other. 



correct, resulting in defects connected by a chain of errors 
- a logical error. 

This allows us to set a maximum edge length that 
is allowed between connections in the minimum weight 
matching algorithm. Instead of creating a completely 
connected graph structure from the classical data, we in- 
stead create multiple smaller subgraphs, with each sub- 
graph having no connections of weight greater than a 
maximum edge parameter mg. As the separation and 
circumference of defects within the lattice determines the 
effective distance of the quantum code rrie can safely be 
chosen such that rrie = ld/2\. This approximation en- 
sures that all error chains throughout the lattice with a 
weight < [d/2\ are connected within the classical pro- 
cessing layer. This classical approximation is defined to 
fail if a single length > [d/2\ error chain occurs during 
computation. This is much less likely than failure of the 
code itself which can occur as the result of numerous, not 
necessarily connected, arrangements of [{d — l)/2\ errors. 

It should be noted that this approximation does not 
neglect all error chains longer than rrie. However, the 
classical error correction data has no knowledge regarding 
the actual path an error chain has taken through the 
lattice. The mg approximation neglects all error chains 
with endpoints separated by > rrie cells. 

By making this approximation, we speed up the clas- 
sical processing for minimum weight as, on average, the 
algorithm will be run on a very sparse graph structure. 



B. Classical simulations as a function of total 
cluster volume and code-distance. 

Given the above classical approximation, we have 
benchmarked the Blossom V algorithm [38 as a function 



10 



of the total volume of the cluster for various distances of 
the quantum code, d and hence mg. 

In order to investigate the classical processing require- 
ments of the system, we will assume a physical error rate 
of the quantum computer. As with previous investiga- 
tions into this system we will be assuming, throughout 
this discussion, that the quantum computer operates at 
a physical error rate of p = 10"^ [23l |25] and the fault- 
tolerant threshold is pth ^ 0.61% [20]. Given this base 



The me approximation to the subroutine prohibits the 
establishment of graph connections between two cluster 
cell co-ordinates separated by a distance > mg = d/2. 



assumption. Figs. 11 and 12 examine the processing time 
of the modified Blossom V algorithm run on a MacBook 
Pro (technical details of this computer are summarized 
in Appendix. [a|. For various volume sizes, from 35^ 
to 295^ cells, random single qubit Z errors were gener- 
ated with a probability of 10~^. The processing time 
for each value of V was examined for different distance 
quantum codes and hence different values of the approx- 
imation parameter, rrie = d/2. A list of cell co-ordinate 
tuples, (i^j^T)^ corresponding to the endpoints of error 
chains (cells with changed parity) was constructed. This 
list is the input to the classical processor as, in practice, 
it would be provided directly by the hardware. 

The tuple data was then processed in two stages. In 
the first stage, which we denote as the tree creation layer, 
an 8-way sort /search tree was created from the tuple in- 
formation and then used to generate a list of connec- 
tions (edges) between cells and their distances (weights) 
[Fig. [ll]. The second stage of the simulation applied the 
Blossom V minimum weight matching algorithm to the 
generated sparse graph structures [Fig. [12] . The bench- 
marking data was taken using 10^ runs per data point. 

For each cluster volume, code distances d = 
[8, 10, 12, 14, 16, 18, 20], corresponding to rrie = [4, .., 10] 
were simulated. For both the tuple tree creation and 
the minimum weight matching. Figs. [TT] and [12] illustrate 
that the approximation parameter, mg, does not signif- 
icantly alter the total simulation time and that for a 
given (y,me), tree creation and minimum weight match- 
ing take a similar amount of time. 



C. Parallelizing the algorithm 



The numerical simulations shown in Figs. 11 and 12 
clearly illustrate that the minimum weight matching sub- 
routine cannot be run over the entire lattice used for 
TCQC. As a rough estimate, a mainframe device such as 
the one introduced in Ref. [25] consists of a lattice cross 
section measuring (5 x 10^) x (4 x 10^) unit cells. Clearly 
in order to achieve classical processing speeds of the or- 
der of microseconds (for any distance topological code), 
either the classical fabrication of the processing equip- 
ment must allow for at least a 10-15 order of magnitude 
speed up from a standard Laptop or the application of 
the tree creation and minimum weight subroutines must 
be highly parallelized. 

Due to the approximation made to the Blossom V al- 
gorithm, parallelizing the classical processing is possible. 



In Fig. 13 we illustrate the relative frequency of differ- 
ent sized connected components within the lattice at an 
error probability of p = 10~^, for rrie = [4, .., 10]. These 
simulations were performed using the Floyd- Warshall al- 
gorithm [43, 44 obtained for a volume region ofV = 50^, 
with 2 X 10^ statistical runs (resulting in approximately 
3 X 10^ connected components in total). In these simu- 
lations, the size of each connected component in the lat- 
tice, n(me), does not represent the longest path through 
the graph. Instead it represents the physical edge length 
through the cluster of a cube of sufficient size to fully 
contain each connected component. In Appendix [B] we 
provide additional simulations showing the distribution 
of connected components within the cluster to assist in 



the explanation of Fig. 13 



Parallelizing the minimum- weight matching procedure 
requires subdividing up a large volume of the cluster into 
smaller regions such that each instance of the tuple 
tree creation and the minimum weight algorithms faith- 
fully return the same results as processing the entire vol- 
ume (up to the failure probability of the computer). In 
Fig. [13] we provide an approximate scaling of the decay 
of each curve, representing the volume independent rela- 
tive frequency of connected component with a linear size, 
n(me), in the cluster. In order to parallelize classical pro- 



cessing, two regions are defined. Fig. 14 i. illustrates. 

The inner volume defines the minimum weight process- 
ing region while the outer volume, with an edge length of 
3x the inner volume, defines the tree creation processing 
region. During tuple tree creation, if any connected 
component contains at least one vertex within the inner 
volume it is sent to independent instances of minimum 
weight matching. Provided that the edge length of these 
regions are large enough, then any connected component 
with at least one vertex in the inner volume will be fully 
contained within the outer volume with high probability. 

To determine the size of these processing regions we 
utilize the decay of the curves in Fig. 13 The prob- 



ability of failure when parallelizing classical processing 
should be approximately the same as the failure proba- 
bility of the quantum computer itself. In order to deter- 
mine these failure probabilities, we consider the volume 
of the cluster required to perform a logical CNOT op- 
eration as a function of rrie. Fig. [Toj illustrates the log- 
ical structure of the lattice. Each logical qubit cell in 
the cluster consists of a cluster cross-section measuring 
{2d + d/2){d + d/4) = 2bml/2 cells. A CNOT gate re- 
quires 4 logical cells and the depth through the cluster 
required to perform the gate is {2d-\-d/2) = brrie cells [20] . 
Hence the total cluster volume for a CNOT operation is 
V 250m^ cells. 

The failure probability of the quantum computer dur- 
ing a logical CNOT is approximately. 
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FIG. 11: Benchmarking data for the tuple tree creation processing layer. Taken with 10^ statistical runs, an operational 
error rate of p = lO"'^ for various sized codes, d = (8, 10, 12, 14, 16, 18, 20), rUe = [4, .., 10]. Notice that rUe does not significantly 
alter the processing time. 



where, 



\pthj 



10 



-2me 



(8) 



is the probability of failure for a single logical qubit a 
single layer thick, with a fault-tolerant threshold of pth ^ 
0.61%, p = 10-"^ and A(me) = 4.{2d + d/2) = 20me is 
the number of such layers of the cluster that need to be 
consumed to perform a logical CNOT operation. 

Given the failure rate of the quantum computer, we 



utilize the data from Fig.[l3]to determine the edge length 
of a volume large enough to encapsulate all connected 
components with a probability approximately equal to 
Eq.[7| As Fig. 13 represent relative frequencies (the prob- 
ability of a connected component of size n(me), relative 
to a connected component of size one) we scale P(n, mg) 
by the number of isolated errors expected in a cluster 
volume required for a CNOT. Hence, 

6p X 250m^P(n, me) = 1 - (1 - l^(me))^^'^^^ (9) 

where the factor of 6 comes from the 6 independent qubits 
per unit cell of the cluster. Eq.|9]is then solved for n(me) 



givmg, 
n(me) 



(3{me) 



1 _ /I _ 1 n — 2me\20me 



where [a{me)^ /3{me)] are the scaling parameters shown 
in Fig. [I3l 

The values of nirrie) for rrie = [4, ..10] and the proba- 
bilities of the logical CNOT failure and equally the prob- 
ability of a connected component of size greater than n 
are shown in Tab. [H 

In Tab.|l]we give the size of the processing regions for 
parallelizing both the tuple tree creation and mini- 
mum weight matching processes such that the probabil- 
ity of any connected component within the inner region 
unbounded by the boundary of the outer region is ap- 
proximately the same as the failure rate of the quantum 
computer. The value n(me) therefore sets the edge length 
of the inner and outer volume regions. 

Using this estimate, the tree creation layer becomes 
an interlaced network, with each individual instance of 
tree creation operating over a volume oiV ^ 27n(me)^. 
As we can neglect larger connected components (which 
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FIG. 12: Benchmarking data for the minimum weight matching processing layer. The simulation conditions for this data set 
are identical to the tree creation layer shown in Fig. |11[ 



rUe 


CNOT failure 


n(me) 


P(connected component > n) 


4 


0(10-^) 


10 


0(10-^) 


5 


0(10-^) 


15 


0(10-^) 


6 


o(io-^°) 


23 


o(io-^°) 


7 


0(10-^') 


32 


0(10-^') 


8 


0(10-^^) 


44 


0(10-^^) 


9 


0(10-^^) 


62 


0(10-^^) 


10 


0(10-^^) 


81 


0(10-^^) 



TABLE I: Maximum edge length, n(me), of a cube of suffi- 
cient size in the lattice to encapsulate all connected compo- 
nents of the tuple tree creation graph structure. The last 
column is the probability that a connected component of the 
graph is unbounded by a cube of volume and Eq.jo] ensures 
that this occurs with approximately the same probability as 
the CNOT failure rate of the topological computer. 



occur with probability roughly equal to the probabil- 
ity of quantum failure), any connected component that 
contains at least one vertex within a central volume 



V ~ n{meY will be fully contained within the outer vol- 
ume oi V ^ 27n(me)^. The central volume represents 
the region that is sent to separate instances of the mini- 
mum weight matching algorithm. The tree creation layer 
is interlaced such that each of the central volumes touch 
but do not overlap. 

This processing structure ensures that multiple paral- 
lel instances of minimum weight matching will produce 
identical results to an individual instance of minimum 
weight matching run over the entire volume. As the tree 
creation layers are interlaced, tree structures that cross 
the boundaries of two inner volumes regions will be sent 
to two independent instances of minimum weight match- 
ing. After processing, these duplicate results will simply 
be removed from the final error list. 

We can now combine the results from Tab. [J with the 
simulation data of Figs . pT] and p!2] to determine the maxi- 
mum size and speed of a quantum computer such that er- 
ror correction data can be processed sufficiently quickly. 
Fig.p!4)3. illustrates how error correction data is collected 
as qubits along the third axis of the cluster are sequen- 
tially measured. As the tuple tree creation processing 
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FIG. 13: Volume independent distribution of connected component sizes for cluster error data. Shown above is the distribution 



of the maximum linear size of each connected component for error data, simulated with p = 10 



= [4, .., 10] and performed 



with 2 X 10^ statistical runs (giving a total number of connected components ^ 3 x 10^). Simulations were performed using 
the Floyd-Warshall algorithm 43 , 44 , but instead of calculating the maximum distance between any two nodes in a connected 
graph we instead calculate the maximum physical distance through the cluster (in a single spatial dimension) between any two 
connected nodes. These results allow us to estimate the edge length of a cube of sufficient size to encapsulate all connected 
components at a given value of rrie. Performing an approximate exponential fit to the decay of these curves allow for the 
estimation of the probability of obtaining very large connected components. Appendix ^ presents further simulation results 
explaining the general properties of this curve. 



layer consists of an interlaced set of V = 27n(me)^ cubes, 
two thirds of the data for any given volume is taken from 
the previously collected results while the final one third 
consists of newly collected data. Therefore the processing 
"window" available for each instance of the tuple tree 
creation (and hence each instance of minimum weight 
matching associated with each tree creation volume) is 
the time required to collect this new data. 



The optical network illustrated in Sec. Ill has each par- 
ity calculation from the detector banks, for each unit cell, 
occurring over three successive cluster faces and, in the 
absence of loss, occur every 2T, where T is the separation 
period of photons (the quantum clock rate). Tree cre- 
ation from the parity tuples for a given volume element 
utilizes the last 9n(me)^ x 2n(me) tuple information from 
the previous instance of tree creation and must store the 
same amount for use in the next tree creation subrou- 



tine. Parity tuple tree creation processing is repeated 
for each 27n(me)^ volume every 2Tn{me) seconds. 
Taking the simulation data from Fig. [iT] (as tuple 
tree-creation is the slower of the two processes), the 
fastest clock cycle, Tmin, can be calculated for each value 
of rrip as, 



i(^e) = 



t(3n(me)) 
2n(me) 



(11) 



where t{3n{me)) is the processing time as a function of 
the edge length of the processing volume, 3n(me), shown 
in Fig. 12 The results are shown in Tab. Ill) 



The last column in Tab. Hll calculates the total number 
of processing instances required per logical qubit in the 
computer. This is calculated as the ratio of the cross- 
sectional area of a logical qubit to the cross-sectional area 
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FIG. 14: Tree creation structure and minimum weight matching processing structure for parallel application of minimum weight 
matching, a) illustrates the volume regions for the two classical processes. For a given approximation parameter, me, a volume 
of 27n(me)^ of tuple data is sent to an independent tree creation process. Once processed, any subgraph with at least one vertex 
in the inner volume region of n{me)^ is sent to a separate instance of minimum weight processing for this region. Parallelization 
is achieved by interlacing the outer volumes such that the inner volumes touch, b) illustrates the temporal processing of this 
data, where a given tree creation volume is processed through time. For each instance of the tree creation process, 67% of the 
tuple data is taken from the previous instance of tree creation (as the previous volume overlaps the new volume by two-thirds) . 
The remaining 33% of the data is collected directly from the measurement of photons. Therefore both the tree creation process 
and minimum weight process must be completed within the time for collection of the new data in order to keep up with the 
quantum clock speed. 



me 


CNOT failure 


n{me) 


Tmin(me)(/is) 


"window", 2n(me)(/is) 


CNOT operating Freq. 


Processing instances / Logical qubit, (I) 


4 


0(10-^) 


10 


0.06 


19 


105 kHz 


8.7 


5 


0(10-^) 


15 


0.28 


31 


18 kHz 


5.3 


6 


o(io-^°) 


23 


1.1 


46 


4 kHz 


3.4 


7 


0(10-^') 


32 


3.2 


64 


1 kHz 


2.4 


8 


0(10-^^) 


44 


9.3 


87 


0.3 kHz 


1.7 


9 


0(10-^^) 


62 


28 


124 


0.1 kHz 


1.0 


10 


0(10-^^) 


81 


68 


162 


37 Hz 


0.77 



TABLE H: Maximum size and speeds for topological quantum computers when classical processing is performed utilizing the 
benchmarking data of Figs. 11 and |12[ The failure of logical CNOT gates defines the size of the computer, ^ 1/KQ, where 
Q is the number of logical qubit s in the system and K is the total number of logical time steps available for an algorithm. 
Tmin{me) dcfines the maximum speed the quantum network can be operated such that error correction data can be processed 
sufficiently quickly. The processing "window", independent of the benchmarking data, is related to the parallelization of 
classical processing. Processing instances / Logical qubit defines how many classical processes are required for a lattice cross 
section housing one logical qubit. 



of the minimum weight matching processing volume. 

/ = 4x(^^^±iP^=4x44, (12) 

n{me) 4n(me) 

the factor of 4 is introduced since each instance of min- 
imum weight matching has an associated tuple tree 
creation process and that primal and dual lattice error 
correction is performed independently giving another fac- 



tor of two. As the size of the quantum computer in- 
creases (increasing mg) this ratio decreases as the scaling 
of nirrie) is approximately n{me) ^ 0{ml). 

While the total size and speed of a topological quan- 
tum computer will ultimately be governed by the exper- 
imental accuracy in constructing each quantum compo- 
nent, the results shown in Tab. [ll] are promising. As- 
suming that quantum fabrication can reach and accu- 



15 



racy of p = 10~^, current classical technology is quite 
sufficient to process error correction data for a large 
range of computer sizes. The logical failure rate of the 
CNOT gate approximately defines the size of the com- 
puter, j9l(CN0T) ^ 0{1/KQ), where Q is the number 
of logical qubits in the computer and K is the number of 
logical time steps in a desired quantum algorithm (note 
that the application of non-clifford gates will lower this 
effective size further). Even for a small topological com- 
puter {rrie = 4 has sufficient protection for 1000 logical 
qubits running an algorithm requiring approximately 10^ 
time steps) less than ten classical processing instances 
are required per logical qubit, a quantum network run at 
^ 17MHz with a logical CNOT operating frequency of 
^ lOOkHz. 

The classical processing power utilized in this inves- 
tigation is clearly not specially designed for the task of 
operating a topological computer. Not only can we safely 
assume that classical processing power will increase be- 
fore the advent of a large topological computer, but the 
design and implementation of both specialized hardware 
and more optimal coding should also result in significant 
increases in the achievable operational frequency of the 
quantum network and logical gates. More recent analysis 
has suggested that possible operational frequency of the 
quantum network could reach the lOOMHz level [45 . In 
this case, if classical processing could result in a 2-3 order 
of magnitude speed up when moving to optimized hard- 
ware and software, current classical technology would be 
sufficient for a quantum computer capable of a logical 
error rate ~ O(10~^^) and a logical CNOT frequency of 
^ 0.3MHz. 



VI. OBSERVATION AND CONCLUSIONS 

In this work we have focused exclusively on the classi- 
cal requirements to perform the underlying error correc- 
tion processing for the network. As the error correction 
procedures can be thought of as the base-level process- 
ing, the development of an appropriate quantum con- 
troller system is an obvious next step. This higher level 
classical controller will essentially be responsible for the 
following, 

1. The compilation of a user-designed quantum cir- 
cuit into an appropriate measurement sequence on 
a topological cluster. 

2. Direct control of quantum components within the 
measurement system of the topological cluster in 
order to change the measurement basis for the pho- 
ton stream. 

3. The dynamic allocation of cluster resources depen- 
dent on operational error rate. Specifically, the 
fundamental partitioning of the lattice into appro- 
priately separated defect regions for logical qubit 
storage. 



4. Accepting the data from the error correction pro- 
cessing layer to faithfully ensure accurate error cor- 
rection is performed during computation. 



5. Dynamical restructuring of the topological lattice 
partitioning to allow for ancilla preparation for non- 
Clifford quantum gates, and optimization of logical 
qubit /qubit interactions for specific quantum sub- 
routines. 



This last point is one of the more interesting questions 
that can be addressed in this model. As we noted in the 
introduction, once the cluster lattice is prepared, data 
processing is performed via software. The structure of 
the topological lattice essentially allows for qubit /qubit 
interactions in a 2D arrangement [Fig. 10 . However, 



provided we have access to a large cluster lattice, we can 
envisage the dynamical creation of data pathways and 
"flying defects" in order to speed up specific quantum 
subroutines. This could lead to some extremely inter- 
esting avenues of investigation in software control and 
optimization of a TCQC architecture. 

This analysis has demonstrated that the classical error 
correction requirements necessary to construct an opti- 
cal quantum computer based on topological cluster states 
is certainly feasible given today's technology. We have 
illustrated how minimum weight matching, required to 
process error data from the topological mainframe, can 
be optimized in such a way to allow for a massively par- 
allelized processing network that can process information 
for large topological clusters. 

These results are very encouraging. As with the quan- 
tum preparation network, the classical front end can also 
be constructed in a modular manner. As the quantum 
preparation network is increased in sized via the addition 
of more photonic chips, the classical processing network is 
also expanded in a similar way. Parity check processors, 
tuple tree creation processors and minimum weight 
matching processors are also linked into a pre-existing 
network as its size expands. The results of this investi- 
gation give us a very optimistic outlook on the viability 
of the topological model as a possible avenue to achieve 
truly large scale quantum computation. 
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APPENDIX A: TECHNICAL SPECIFICATIONS 
FOR SIMULATIONS. 



APPENDIX B: SIMULATIONS OF GRAPH SIZES 
FOR CLUSTER ERRORS. 



The technical specifications to the computer used in 
benchmarking simulations are summarized below. 



Process 


Benchmark 


Floating Point Basic 


3.1Gflop/s 


vecLib FFT 


3.5 Gflop/s 


Memory Fill 


6.2 GB/s 



TABLE IIL Benchmarking for system processes for the Mac- 
Book Pro 2.2 GHz, 3GB RAM Benchmarking data was taken 
using the program XBench, version 1.3. 



The simulations shown in Fig. 13 illustrate the distri- 
bution of the largest physical distance through the clus- 
ter (in a single spatial dimension) between any two nodes 
for each connected component for graph structures estab- 
lished using the rrie approximation detailed in the main 
text. The following results illustrate the structure of this 
distribution in more detail. 
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FIG. 15: The maximum graph diameter over 10^ statistical runs utihzing p = 10~^, rrie = 6 and a total cluster volume of 
V = 100^. The Floyd- Warshall algorithm is designed to find the shortest pathway between any two connected nodes in the 
cluster and in these simulations we maximize over all possible connections. Unlike Fig. ^] we are examining the actual graph 
diameter, D, and not the linear separation of nodes in the physical cluster. Here we are simply finding the maximum graph 
diameter in the complete data set, unlike Fig. 13 which calculates the diameter of all connected components (hence these results 
exhibit volume dependence). At p = 10~^ and nie = 6, the graph structure for D = 8 is the most probable. Changing p shifts 
which diameter graph is the most probable, while changing nie changes the values of D where peaks occur. 
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FIG. 16: Volume independence of distribution of all connected components within a cluster volume. The upper four curves 
are simulations calculating the largest connected component in cluster volumes of V = (50^,75^, 100^, 125^) while the lower 
four curves examine the distribution for all connected components at p = lO"'^ and irie = 6 (total number of simulations 
vary between O(IO^) - O(IO^)). As you can see, volume independence when calculating all connected components is good. 
Additionally, we are now calculating the maximum p%5zca/ separation of endpoints within each connected component. This has 
the effect of smoothing out the curves in Fig. 15 and slightly shifting the main peak to the left (as maximum graph diameter, 
in general, is larger than the maximum node separation in the physical cluster). In the simulations calculating the largest 
connected component, the main peaks shift to the right as volume increases. This is again due to the fact that the largest 
connected component will scale with volume. 



